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ABSTRACT : Genetic diversities, population genetic structures and demographic histories of the thread-sail filefish Stephanolepis 
cirrhifer were investigated by nucleotide sequencing of 336 base pairs of the mitochondrial DNA (mtDNA) control region in 111 
individuals collected from six populations in Korean coastal waters. A total of 70 haplotypes were defined by 58 variable nucleotide 
sites. The neighbor-joining tree of the 70 haplotypes was shallow and did not provide evidence of geographical associations. Expansion 
of S. cirrhifer populations began approximate 51,000 to 102,000 years before present, correlating with the period of sea level rise since 
the late Pleistocene glacial maximum. High levels of haplotype diversities (0.974±0.029 to 1 .00010.076) and nucleotide diversities 
(0.014 to 0.019), and low levels of genetic differentiation among populations inferred from pairwise population F sx values (-0.007 to 
0.107), support an expansion of the 5. cirrhifer population. Hierarchical analysis of molecular variance (AMOVA) revealed weak but 
significant genetic structures among three groups (F cx = 0.028, p<0.05), and no genetic variation within groups (0.53%; F sc = 0.005, p 
= 0.23). These results may help establish appropriate fishery management strategies for stocks of S. cirrhifer and related species. 
(Key Words : Stephanolepis cirrhifer, Mitochondrial Control Region, Demographic History, Population Genetic Structure, Nucleotide 
Sequence Analysis) 



INTRODUCTION 

Filefish (family Monacanthidae) include 95 species 
widely distributed in both temperate and tropical seas 
(Nelson, 1994). The thread-sail filefish, Stephanolepis 
cirrhifer (Temminck and Schlegel), is an economically 
important species for the Korean fishery and aquaculture 
industries. S. cirrhifer aquaculture has received particular 
attention because of the fish's high per weight market price 
and rapid growth rates (market size is reached in only one 
year) (Miyajima et al., 2011). Although, the species is also 
widely distributed in the western Pacific Ocean, occurring 
in the East Sea, the Yellow Sea, and the East China Sea 
(Masuda et al., 1984; Shao et al., 1990; Nelson, 1994; Ni 
and Kwok, 1999; Kim et al., 2005) most commercial 
catches are from the southern coastal area of the Korean 
Peninsula and the Kuroshio Current of the tropical Pacific 
Ocean. 

In 1985, the S. cirrhifer catch in the Northwest Pacific 
exceeded 250,000 tonnes, but it had decreased to less than 
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950 tonnes by 2002 (FAO, 2004). Stocks of filefish species 
are decreasing as a result of overfishing and environmental 
degradation. Therefore, stock enhancement programs are 
critically needed to increase severely depleted fishery 
stocks. Hatchery-based release of filefish is one approach to 
revitalizing stocks and yields. However, hatchery stocks are 
susceptible to genetic degradation and reduction of genetic 
diversity: release of hatchery stocks can threaten the 
integrity of wild populations by changing the genetic 
composition and diversity of the gene pool (Allendorf and 
Ryman, 1987). Thus, evaluation of the biological and 
genetic characteristics of natural S. cirrhifer populations 
requires urgent attention to maintain genetic diversity levels 
and promote sustainable harvests of natural populations. 

The primary aim of many historical and evolutionary 
studies of marine species is to ascertain oceanographic 
conditions, climatic factors, and vicariant events disrupting 
gene flow and genetic diversity, in relation to biogeographic 
patterns of demographic expansion and contraction. Genetic 
diversity is important for the long-term persistence and 
survival of natural populations because it confers the ability 
to adapt to changes in environmental conditions (Frankel 
and Soule, 1981). The spatial distributions of population 
genetic structures are influenced by the effects of gene flow, 
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natural selection, genetic drift and mutation (Slatkin, 1987), 
which may be affected by both current environmental 
conditions (Weersing and Toonen, 2009) and historical 
biogeographical processes (Avise, 2000). Assessments of 
the genetic structure of fish populations can provide 
invaluable information for successful conservation or 
management of exploited species (Bailey, 1997; Zang et al., 
2006). 

Information on nucleotide sequences in maternally 
inherited mitochondrial DNA (mtDNA) genes is useful for 
investigating interspecific and intraspecific genetic diversity 
among closely related taxa, species, or populations (Wilson 
et al., 1985; Moritz et al., 1987; Avise, 1994). The 
intraspecific population genetic structures of S. cirrhifer 
population are still poorly understood. A recent study by 
Yoon et al. (2011) found no differences in the nucleotide 
sequences of the mtDNA cytochrome b gene (Cyt b) 
between S. cirrhifer populations off the southern and 
eastern coasts of the Korean Peninsula. They suggest that 
substantial gene flow occurs within and between 
populations of S. cirrhifer, and that western Pacific S. 
cirrhifer belong to a single panmictic population; however, 
this study has not investigated the genetic structure of 
populations from the west coast of Korea. 

In the mtDNA genome, the sequence variability of the 
non-coding control region is often higher than the 
variability in coding regions; thus, sequencing of the non- 
coding region has been recommended for investigating 
genetic diversity at the intraspecific level (Moritz et al., 
1987). Up until now, the mtDNA control region has been 
generally used to study genetic diversities, population 
structures, and intraspecific phylogenesis of fish (Brown et 
al., 1993; Stepien and Faber, 1998; Sato et al., 2001; 
Guarnieo et al., 2002). 

In the present study, we examine the genetic diversity 
and population structure of S. cirrhifer along the Korean 



coast using sequence analysis of the 5' end of the mtDNA 
control region to obtain basic data for fisheries resources 
management. 

MATERIALS AND METHODS 

Sampling and sequencing 

Stephanolepis cirrhifer specimens were collected in 
2009 and 2010 from six sites located on three Korean seas: 
the Yellow Sea (one site), South Sea of Korea (South Sea) 
(three sites), and East Sea (two sites) (Table 1). All the 
samples examined, except the Yellow Sea, for the present 
mtDNA (control region gene) analysis had been used 
previously for the mtDNA analysis (Yoon et al., 2011). 
Specimens were stored at -80°C or preserved in ethanol at 
room temperature until DNA extraction. DNA was extracted 
from approximately 70 mg of each specimen with a 
QIAGEN Blood and Cell Culture DNA Midi Kit (Qiagen, 
Germany) according to the manufacturer's instructions. 

The polymerase chain reaction (PCR) was used to 
amplify the mtDNA control region with the primers 
StepCR-F (5 ' -CTAGCTCCC A AGCTAGG ATT-3 ') and 
StepCR-R (5'-TGGTGAGCCACGTATTGCAA-3'), newly 
designed on the basis of the complete mitogenomic 
sequence of S. cirrhifer, available in the GenBank genetic 
sequence data bank (NC_003177). PCR amplification was 
performed with a DNA Engine thermocycler (MJ Research, 
Tokyo, Japan) in 20-|il reaction tubes containing 100 ng of 
genomic DNA, 2 |iM of each primer, 0.25 mM of each 
dNTP, 1 unit of TaKaRa LA Taq DNA polymerase (Takara 
Shuzo, Shiga, Japan), and 2 of 10x LA Taq reaction 
buffer (Takara Shuzo). The PCR conditions consisted of 
preheating at 94°C for 5 min, followed by 35 cycles of 
94°Cfor 30 s, 58°Cfor 30 s, and 72°C for 30 s, with a final 
7 min extension at 72°C. 



Table 1. Sampling sites, dates, geographical coordinates, number of individuals examined, and measures of mtDNA diversity in six 
Stephanolepis cirrhifer populations 



Sampling site 


Abbreviation 


Date of 


Geographical co-ordinates 


N 


No. of 


h 


7t 


collection 


Latitude 


Longitude 


haplotypes 


Yellow Sea 


















Wi-do 


WID 


Oct. 2010 


35°35'33.14"N 


126°16'23.65"E 


1 


7 


1.000±0.076 


0.018 


South Sea 


















Jeju 


JEJ 


Aug. 2010 


33°03'14.54"N 


126°32'39.67"E 


18 


15 


0.974±0.029 


0.014 


Yeosu 


YES 


Aug. 2010 


34°41'22.99"N 


127°39'59.62"E 


23 


21 


0.992±0.015 


0.014 


Geoje 


GUJ 


July. 2010 


34°43'15.11"N 


128°41'07.46"E 


24 


18 


0.967±0.024 


0.015 


East Sea 


















Pohang 


POH 


Oct. 2010 


36°02'50.30"N 


129°24'31.75"E 


18 


15 


0.980±0.024 


0.014 


Uljin 


ULJ 


Oct. 2009 


36°58'04.20"N 


129°26'09.39"E 


21 


18 


0.986±0.019 


0.019 



N: Sample size, h: Haplotype diversity, %: Nucleotide diversity. 
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The size of the PCR product was verified by 1.0% 
agarose gel electrophoresis after ethidium bromide staining. 
The PCR product was purified with the AccuPrep PCR 
Purification Kit (Bioneer, Daejeon, Korea). After cycle 
sequencing with the ABI PRISM BigDye Terminator v3.1 
Cycle Sequencing Ready Reaction Kit (Applied Biosystems, 
Foster City, CA, USA), the purified PCR product was 
sequenced directly on an ABI 3730x1 DNA Analyzer 
(Applied Biosystems). The primers used for sequencing 
were the same as those used for PCR amplification. 



determined. The evolutionary rate for the control region 
seems to vary among major taxonomic groups of marine 
fishes. For the Japanese sea bass {Lateolabrax japonicas) 
and spotted sea bass (Lateolabrax maculates), the sequence 
divergence rate of the control region is 6% per million years 
(Liu et al., 2006). In Nibea albiflora, the sequence 
divergence rate in the control region is 5% to 10% per 
million years (Han et al., 2008). In the present study, we 
calculated a sequence divergence rate of 5% to 10% per 
million years for the control region of S. cirrhifer. 



Data analysis 

The sequence fragments obtained in this study were 
aligned with GENETIX-WIN ver. 4.0.1 (Software 
Development Co., Ltd, Tokyo, Japan) to identify sequence 
variants. The integrated software package DnaSP v. 4.90.1 
(Rozas and Rozas, 1997) was used to determine haplotypes. 
Haplotype diversity (h) and nucleotide diversity (n) were 
estimated according to Nei (1987), based on Kimura's two- 
parameter distance method, using the K and DA in the 
REAP program (McElroy et al., 1993). Genetic 
relationships among haplotypes were reconstructed using 
the neighbor-joining method, generated with the Seqboot, 
Neighbor, and Consensus options in PHYLIP v. 3.6 
(Felsenstein, 1993). A bootstrap analysis of 1,000 replicates 
evaluated support for phylogenetic relationships after 
construction of a genetic distance matrix based on 
nucleotide divergences between haplotypes, estimated 
according to Nei (1987) and the Kimura two-parameter 
model (Kimura, 1980). 

Pairwise population F ST values were calculated to 
estimate genetic differentiation between populations 
(according to Slatkin, 1991) using ARLEQUIN v. 3.1 
(Excoffier et al., 2005). The statistical significance of F ST 
values was tested using 10,000 random permutations. The 
variance component was evaluated at each hierarchical 
level using analysis of molecular variance (AMOVA) 
(Excoffier et al., 2005) to assess genetic divergences within 
and among groups, using ARLEQUIN v. 3.1. 

Neutral expectations and historic demographic 
expansions were investigated by examining Fu's F s and 
Tajima's D mismatch distributions with the sudden 
expansion model (Rogers and Harpending, 1992). A 
goodness-of-fit test was used to test the validity of the 
sudden expansion model using a parametric bootstrap 
approach based on the sum of squared deviations (SSD) to 
compare the observed and the estimated mismatch 
distributions (Schneider and Excoffier, 1999). Both the 
neutrality test and the mismatch distribution analysis were 
performed in ARLEQUIN v. 3. 1 (Excoffier et al., 2005). 

The mutation rate of the S. cirrhifer control region gene 
over the estimated time since expansion has not been 



RESULTS AND DISCUSSION 

Genetic diversity 

Analysis of a 336 base pair fragment of the 5'- end of 
the mtDNA control region in 104 individuals of S. cirrhifer 
from six populations defined a total of 70 haplotypes. The 
neighbor-joining tree of the 70 haplotypes was shallow and 
did not provide evidence for geographical associations 
(Figure 1). This suggests a signature of population 
expansion and/or high gene flow among populations. 

Haplotype diversities (h) and nucleotide diversities (n) 
were high for all populations, ranging from 0.974+0.029 
(JEJ) to 1.000±0.076 (WID), and 0.014 (JEJ, YES, POH) to 
0.019 (ULJ), respectively (Table 1). Grant and Bowen 
(1998) suggest four basic scenarios for population 
demographic histories of marine fishes with different 
measures of haplotype and nucleotide diversities. They 
suggest that high haplotype and nucleotide diversities may 
be attributed to a long evolutionary history in a large 
population; examples include the Japanese anchovy (h = 
0.91, n = 0.01) and the Atlantic bluefish (h = 0.70, n = 0.01). 
This condition might also be observed in re-mixed 
populations of individuals from historically separated 
populations (Avise, 2000). Therefore, genetic diversity 
values of S. cirrhifer populations might reflect strong 
dispersal capacities among ancestral populations with a 
long evolutionary history in Korean waters. 

High genetic diversities within and between populations 
provide a potential genetic resource for future adaptation 
(Savolainen and Kuittinen, 2000; Hurt and Phillip, 2004). 
Recent research has examined genetic diversity in other 
fishes to determine the present status of populations in this 
area (e.g., Liu et al., 2006; Shui et al., 2009; Xiao et al., 
2009). 

Population genetic structure 

Estimation of population structures using molecular 
markers is a powerful approach to understanding the 
dynamics of natural populations. Pairwise population F ST 
estimates are generally low (Table 2), suggesting little 
genetic differentiation between population pairs, perhaps 
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Figure 1. Neighbor-joining tree for 70 control region haplotypes of Stephanolepis cirrhifer. Bootstrap supports of >50% in 1,000 
replicates are shown. 



Table 2. Pairwise F ST (below diagonal) and p (above diagonal) values among populations of Stephanolepis cirrhifer 


Population 


WID 


JEJ 


YES 


GUJ 


POH 


ULJ 


WID 




0.05 


0.01 


0.02 


0.01 


0.01 


JEJ 


0.051 




0.78 


0.13 


0.02 


0.08 


YES 


0.074 


-0.014 




0.15 


0.08 


0.35 


GUJ 


0.103 


0.020 


0.016 




0.23 


0.56 


POH 


0.107 


0.043 


0.027 


0.011 




0.38 


ULJ 


0.089 


0.025 


0.004 


-0.007 


0.002 
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Table 3. Results of the hierarchical analysis of molecular variance (AMOVA) based on mtDNA control region sequence data for 
Stephanolepis cirrhifer 



Source of 
variation 



df 



Sum of 
squares 



Variance 
components 



Percentage of 
variation 



F index 



Analysis 1 . One gene pool 

Among populations 5 18.321 0.06127 Va 2.34 

Within populations 105 267.931 2.55172 Vb 97.66 

Analysis 2. Three gene pools (Yellow Sea, South Sea, East Sea) 

Among groups 2 9.797 0.07311 Va 2.77 

Within groups 3 8.524 0.01393 Vb 0.53 

Within populations 105 267.931 2.55172 Vc 96.70 



:0.02 



Fct= 0.028 
F sc = 0.005 
F ST = 0.033 



<0.05 



<0.05 
0.23 
<0.05 



attributable to high gene flow. High gene flow is 
characteristic of marine fishes with passively dispersed 
planktonic larvae and adult migration abilities; their 
dispersals may be strongly affected by sea currents, and 
population substructures are limited (Han et al., 2008; Xiao 
et al., 2009; Kim et al., 2010). The Tsushima Warm Current 
(TWC) diverges from the Kuroshio Current, with the main 
flow entering the East Sea west of the Korean Peninsula and 
a subsidiary flow entering the Yellow Sea (Senjyu, 1999; 
Ichikawa and Beardsley, 2002). Thus, the TWC may 
transport larvae into the southern part of the East Sea, near 
Pohang and Uljin on the Korean coast, and to the Yellow 
Sea, near Wi-do. Previous studies based on mtDNA Cyt b 
nucleotide variation have suggested that the South Sea and 
East Sea S. cirrhifer populations comprise a single 
population (Yoon et al., 2011). Our mtDNA control region 
results also showed that estimates of the population 
differentiation seem to provide similar patterns of genetic 
differentiation among the South Sea and East Sea 
populations. However, pairwise population F ST estimates 
between Wi-do and all other populations are relatively large 
(0.051 to 0.107), with significant p values (<0.05), 
compared to values for all other population pairs (-0.007 to 
0.043) (Table 2), suggesting that significant population 
subdivisions exist at small spatial scales. 

The genetic structure of S. cirrhifer populations around 
the Korean Peninsula was estimated by AMOVA (Table 3). 
Overall genetic variation among six populations (Analysis 
1) was 2.34% (F(jy = 0.02, p<0.05), suggesting the 
possibility of substructures among populations. Analysis of 
variation within three groups (Yellow Sea, South Sea, and 
East Sea; Analysis 2) indicated no genetic variation within 
groups (0.53%; F sc = 0.005, p = 0.23), 96.7% of the total 
variation due to differences within populations (F ST = 0.033, 
p<0.05), and 2.77% of the total variation due to differences 
among groups (F CT = 0.028, p<0.05). The hierarchical 
pattern of genetic differentiation among groups of S. 
cirrhifer in Korean waters indicates weak but historical 
patterns of isolation and restriction of gene flow between 
groups. These results are consistent with previous results 



based on the mtDNA Cyt b sequence data (Yoon et al., 
2011), suggesting that the sequence variability of the 
mtDNA control region of S. cirrhifer could be expected to 
provide a powerful means with an increased accuracy and 
resolution to reveal genetic variation as compared with the 
those of the Cyt b gene. However, differences between the 
three groups likely arose during isolation into refugia 
during glacial periods. Rare migrations of strays from the 
various refugia may have occurred during glaciations. 
Pleistocene glaciations (2.5 Ma to 10,000 years ago) 
influenced not only historical demographic patterns (glacial 
population extinctions and interglacial re-colonizations) but 
also the demographics of contemporary populations of 
phylogroups of marine species (Seeb and Crane, 1999; 
Kitano et al., 2007). 

At present, the East Sea is a semi-enclosed marginal sea, 
which was geographically isolated from the Pacific Ocean 
during the Pleistocene glacial period (Nishimura, 1974), by 
a large land bridge that extended from the Yellow Sea to 
Taiwan (Kimura, 1996, 2000). The present population 
structure of S. cirrhifer populations may represent the 
influence of this vicariant barrier. Populations of S. cirrhifer 
may have been isolated in the Yellow Sea, the South Sea, 
and the East Sea during times of glacial maxima. 

Historic demography 

Neutral expectations and historical demographic 
expansions of S. cirrhifer were investigated using Fu's F s 
and Tajima's D tests, and mismatch distributions with the 
sudden expansion model. Fu's F s (-25.427; p = 0.000) and 
Tajima's D (-1.708; p = 0.020) showed significant 
deviations from the neutral evolution model. The mismatch 
distributions of S. cirrhifer in Korean coastal waters 
appeared to be unimodal (Figure 2). Significant differences 
for the sums of the square deviations (p SSD <0.05) between 
the observed and simulated mismatch distributions indicate 
that the population is at equilibrium (i.e., a non-expansion 
phase). The goodness-of-fit test (Pssd = 0.657) did not reject 
the null hypothesis of sudden population expansion for S. 
cirrhifer populations. Thus, both the neutrality tests and the 
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Figure 2. Mismatch distribution constructed using pairwise differences among the mtDNA haplotypes of Stephanolepis cirrhifer. The 
bars are observed pairwise differences values and the open circles and solid line are the expected mismatch distributions under the 
sudden expansion model. 



mismatch distributions indicate a recent expansion of S. 
cirrhifer populations. 

Using x (3.436), the expansion was estimated to have 
occurred 51 to 102 Ma before present, which is consistent 
with sea level rise since the late Pleistocene glacial 
maximum. Rising sea level removed barriers to migration 
and dispersal during the period of postglacial warming, 
establishing potential gene flow between previously 
isolated populations. High levels of genetic diversity and 
low levels of genetic differentiation determined in the 
present study also support population expansion since the 
late Pleistocene. 

In conclusion, this study found high levels of gene flow 
among populations of S. cirrhifer, weak but significant 
genetic differentiation of populations on the west coast of 
Korea (Yellow Sea) and pooled populations in southern and 
eastern Korean coastal waters. Incorporation of additional 
population in the Yellow Sea thus provided an increased 
resolution in the geographic differentiation among three 
groups of populations. Moreover, the genetic structures of S. 
cirrhifer populations indicate the presence of three 
geographic population groups (Yellow Sea, South Sea, and 
East Sea), each of which may require separate fisheries 
management strategies. 
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